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Abstract. We model the cosmological co-evolution of galaxies and their central supermas- 
sive black holes (BHs) within a semi-analytical framevyork developed on the outputs of the 
Millennium Simulation jCroton et al.l2006l : lDe Lucia & Blaizotl2007h . In this work, we an- 
alyze the model BH scaling relations, fundamental plane and mass function, and compare 
them with th e most recent obser vational data. Furthermore, we extend the original code 
developed bv lCroton et alj l l2006h to follow the evolution of the BH mass accretion and its 
conversion into radiation, and compare the derived AGN bolometric luminosity function 
with the observed one. We find, for the most part, a very good agreement between predicted 
and observed BH properties. Moreover, the model is in good agreement with the observed 
AGN number density in < z < 5, provided it is assum ed that the cold gas fraction accreted 
by BHs at high redshifts is larger than at low redshifts dMarulli et alj2003) . 
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1 . The model a set of recipes to follow the BH accretion his- 

tory and the AGN phenomenon. This model 
Our cosmological model for the co-evolution is described in deta il in iCrotoneta l. (2006), 
of DM haloes, galaxies and their central BHs |De Lucia & BlaT^ (l2007h and iMaruUi et all 



consists of three ingredients: i) a numerical N- ^20081). In the f^wing, we just give a brief 



body simulation, the Millennium Run, to ob- description of the assumptions introduced to 
tain the merger history of the DM haloes and describe the BH and AGN evolution, 
subhaloes within the framework of the ACDM 
model, ii) a set of analytic prescriptions to trace 
the formation and evolution of galaxies and iii) 
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1.1. Black Holes and Active Galactic 
Nuclei 

1.1.1. Quasar Mode 

As it is well established that galaxy major 
mergers cannot constitute the only trigger to 
accretion ep isodes in the loca l BH popula- 
tion (see e.g. Marulli et al.''2007', and reference 
therein), we assume that BHs can accrete mass 
after every galaxy mergers, both through co- 
alescence with another BH and by accreting 
cold gas, the latter being th e dominant accre- 
tion mechanism (see e.g. Maru lli et al.l l2006l 
and reference therein). For simplicity, the BH 
coalescence is modelled as a direct sum of the 
progenitor masses, thus ignoring gravitational 
wave losses. We assume that the gas mass ac- 
creted during a merger is proportional to the to- 
tal cold gas mass present, but with an efficiency 
which is lower for smaller mass systems and 
for unequal mergers. This kind of accretion is 
also closely associated with starbursts, which 
occur concurrently. We do not model feedback 
from the quasar activity, but it can be approx- 
imately represented by an enhanced effective 
feedback efficiency for the supemovae associ- 
ated with the intense starburst. 

The evolution of an active BH is described 
as a two-stage process of a rapid, Eddington- 
limited growth up to a peak BH mass, pre- 
ceded and followed by a much longer quiescent 
phase with lower Eddington ratios. In this lat- 
ter phase, the average time spent by AGN per 
logarithmic luminosity interval ca n be approx- 
imated as in i Hopkins et aP (l2005h . 

As discussed in details in iMaruUi et alJ 
(l2008h . with th e original semi-analy tic recipes 
introduced by ICroton et all (|2006|) to follow 
the BH mass accretion, the model underesti- 
mates the number density of luminous AGN 
at high redshifts, independently of the assump- 
tions introduced to describe the accretion effi- 
ciency, the Eddington factor, or the BH seed 
masses. Significant improvement can be ob- 
tained by simply assuming an accretion effi- 
ciency that increases with the redshift. In the 
following Section, we will focus on the predic- 
tions of this new model, which we called 'best 
model' in iMaruUi et all (l2008h . 



1.1.2. Radio Mode 

We assume that, when a static hot halo has 
formed around a galaxy, a fraction of the 
hot gas continuously accretes onto the cen- 
tral BH, causing a low-energy 'radio' activ- 
ity in the galaxy centre. This accretion rate 
is typically orders-of-magnitude below the 
Eddington limit, so that the total mass growth 
of BHs in the radio relative to the quasar mode 
is negligible. It is also assumed that the radio 
mode feedback injects energy efficiently into 
the surrounding medium, which can reduce or 
even stop the cooling flow in the halo cen- 
tres. In this scenario, the effectiveness of radio 
AGN in suppressing cooling flows is greatest 
at late times and for large values of the BH 
mass, which is required to successfully repro- 
duce the luminosities, colours and clustering of 
low-redshift bright galaxies. 

2. Models vs. Observations 

In this Section, we show the most interest- 
ing results about the comparison between our 
model predictions and several observed statis- 
tical properties of BHs and AGN. In all the 
Figures, black and grey symbols show the ob- 
servational data, while red and blue ones show 
the model predictions. 

Figure[T]shows the correlation between the 
masses of the model BHs with six properties 
of their hosts: the K- and B-band bulge mag- 
nitude (Mb and Mk) (Marconi & Hunt 2003i) . 
the bulg e mass and velocity dispersion (Mhuiyp 
and CT r) jHaring & Rixl20"04tlFerrarese & FordI 
l2005h . the circular velocity of the galaxy 
and the virial mass of the DM hal o (Vr 
and Mr.M) dFerraresd l2002t iBaes et all 120031; 
IShankar et al.ll2006h . 

In Figure|2l we compare the BH fundamen- 
tal plane relation in the redshift range 0.1 < z < 
5 predicted by our mod el with that obtained 
bv lHopkins etal] (1200 7a'), using both observa- 
tional data and the outputs of hydrodynamical 
simulations of galaxy merger. 

In Figure [3] we show the BH mass func- 
tion pred icted by our model with those ob- 
served bv'Shanka r et alj ( |2004|) and by Shankar 
(private communication) at z ~ 0. In nei- 
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Fig. 1. Scaling relations between the masses of 
the central BHs in the simulated galaxies with 
six different properties of their hosts: the K- 
and B-band bulge magnitude (top left and right 
panels, respectively), the bulge velocity disper- 
sion and mass (central left and right panels, re- 
spectively), the circular velocity of the galaxy 
(bottom left panel) and the virial mass of the 
DM halo (bottom right panel). Red and black 
dots represent model predictions and observa- 
tions, respectively. Solid blue and dashed black 
and grey lines show the best fit to the model 
predictions and to the observational datasets, 
respectively. 

ther ca se the BH mas ses were determined di- 
rectly: IShankar et aP (2004) derive the BH 
mass from the observed Mbh - i^buige relation, 
while Shankar (private communication) us es 
the Mbh - cr^ relation of iTundo et"aD (l2007h . 

Finally, in Figure |4] we compare the AGN 
bolometric luminosity function predicted by 
our model with several observed ones in differ- 
ent bands. The bolometric corrections adopted 
and the datasets considered ar e the ones dis- 
cussed in lHopkins et al.l (120061) . 

3. Conclusions 

The main results of our analysis are as follows: 
{i) The semi-analytic model is able to re- 
produce the observed BH scaling relations over 
the whole range of BH masses and galaxy 



Fig. 2. The BH fundamental plane in the red- 
shift range 0.1 < z < 5. The red dots are the 
model outputs, while the blue solid lines show 
the best fits to them. The black dashed lines 
show the predictions of Hopkins et al. ( 2007a). 
The galaxy stellar mass, Mjj, is given in units 
of 1O"M0, while the bulge velocity dispersion, 
cr2oo, is in units of 200 km i"' . 



properties probed by observations. The intrin- 
sic scatter in the model is significantly larger 
than in the data, a mismatch that can be ac- 
counted for by adopting the observational se- 
lection criteria to obtain a mock BH catalogue 
with similar characteristics as the observed 
one. 

{ii) We find evidence that a quadratic re- 
lationship provides a significantly better fit to 
some of the model scaling relationships than 
a hne ar one, as already noticed by WvithJ 
(I2006h . 

( Hi) Our model also matches the BH funda- 
mental plane relation derived by lHopkins et al.l 
(l2007ah . and predicts very little evolution of 
this plane, at least out to z ~ 3. 

( iv) The model BH mass function is in good 
agreement with the observed one within the 
mass range accessible by observations, except 
on the range ~ 10^ - lO^M©, in which the num- 
ber density predicted by the model is smaller 
than the observed one. 

(v) Model predictions for the BH mass 
function, scaling relations and fundamental 
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Fig. 3. The model BH mass function (red line) 
co mpared with the one o bservationally derived 
by IShankar et alj (|2004|) (dark grey area), and 
with the new one obtained by Shankar (pri- 
vate communication) (light grey area) using 
the Mbh - o" relation bv lTundo et alj (|2007|) . 



log{koi [^rg sec ']) 

42 44 46 4B 42 44 46 49 42 44 46 48 













; "I 









10 12 14 8 10 12 14 8 10 12 14 

i°e(k.i lk.i.eJ) 



Fig. 4. The model AGN bolometric luminosity 
functions (red lines) compared with several ob- 
served ones (grey dots). The bolometric correc- 
tions adopted and th e datasets considered are 
the ones discussed in lHopkins et alJ (l2006h 



plane relation are basically unaffected when 
using different prescriptions for the AGN 
lightcurves of individual quasar events. 



( vi) The model underestimates the number 
density of luminous AGN at high redshifts, in- 
dependently of the lightcurve model adopted. 
We were not able to eliminate this mismatch by 
simply modifying the accretion efficiency, the 
Eddington factor or the BH seed mass (when 
considered in physically plausible ranges). A 
simple, ad hoc increase of the mass fraction ac- 
creted during the quasar mode at high redshifts 
can indeed remedy the problem. However, this 
solution is not unique as several high-redshift 
modifications to the original model, like new 
mechanisms that trigger BH activity in addi- 
tion to galaxy merging or more efficient gas 
cooling resulting in a larger reservoir of cold 
gas, can be advocated to bring the predictions 
in line with observations. However, it remains 
to be seen whether any of these alternatives is 
physically plausible. 
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